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We have performed total-energy density-functional cal- 
culations using first-principles pseudopotentials to determine 
the atomic and electronic structure of neutral surface and 
subsurface vacancies at the GaP (110) surface. The cation as 
well as the anion surface vacancy show a pronounced inward 
relaxation of the three nearest-neighbor atoms towards the va- 
cancy while the surface point-group symmetry is maintained. 
For both types of vacancies we find a singly occupied level at 
midgap. Subsurface vacancies below the second layer display 
essentially the same properties as bulk defects. Our results for 
vacancies in the second layer show features not observed for 
either surface or bulk vacancies: Large relaxations occur and 
both defects are unstable against the formation of antisite va- 
cancy complexes. Simulating scanning tunneling microscope 
pictures of the different vacancies, we find excellent agreement 
with experimental data for the surface vacancies and predict 
the signatures of subsurface vacancies. 

PACS numbers: 68.35.Dv, 71.15.Mb, 73.20.Hb 



I. INTRODUCTION 

Native defects at and near surfaces are commonly in- 
voked in explaining various electrical and structural prop- 
erties of surfaces and interfaces such as Fermi-level pin- 
ning at clean and metal covered surfaces .Eft] Also prop- 
erties such as the surface chemical reactivity, nucleation 
processes, or diffusion at surfaces and interfaces!] depend 
to a large extent on the presence and properties of native 
defects. Motivation for the study of surface defects also 
arises from the investigation of bulk impurities: Cleav- 
ing crystals in ultrahigh vacuum offers the opportunity 
to study defects at and near the cleavage plane. nPro 
vious investigations on dopant atoms|3~l3 antisitesja and 
defect complexes^ have shown that the concentration of 
such impurities corresponds well to the measured bulk 
concentrations. However, it is by no means obvious how 
the structural and optical properties of the defects are 
changed if they are close to or within the surface. We 
therefore compare theoretically the properties of defects 
at and below surfaces with their corresponding bulk de- 
fects in order to interpret the experimental data mea- 
sured on cleavage planes. 

The aim of the present paper is to analyze how the 
properties and behavior of native defects change with 
their distance from the surface. In particular, we focus 



on vacancies since these .defects have been experimentally 
studiedJja bulk systemsli9~cJ as well as at the (110) sur- 
face. EjU The atomic as well as the electronic structure 
of cation and anion vacancies in different surface layers 
and in a bulk environment are discussed. We also present 
results concerning the formation energies of the different 
vacancies as well as theoretical STM pictures that can be 
directly compared to experiment. We will show that sub- 
surface vacancies are significantly different from surface 
or bulk vacancies: Vacancies in the second surface layer 
are unstable against the formation of antisite vacancy 
pairs. 

The paper is organized as follows: First we briefly de- 
scribe our method. In Sec. Ill we present the atomic and 
electronic structure for gallium and phosphorus surface 
vacancies. In Sec. [V these results are compared with 
calculations of vacancies below the surface and in the 
bulk. We discuss further the formation energies and con- 
clude on the abundance in thermal equilibrium in Sec.^. 
Based on the calculated atomic and electronic structure, 
STM images of the various surface and subsurface defects 
are simulated and will be discussed in Sec. VI. Finally we 
will summarize our results. 



II. CALCULATIONAL METHOD 

Our calculations employ density-functional theory 
using the local-density approximation (LDA) for the 
exchange-correlation functional E»Hp Fully separable, 
norm-conserving pseudopotentialsEJ are used to describe 
the interaction of electrons and ions. The eigenfunc- 
tions p£ the Kohn-Sham operator are expanded in plane 
wavesEll using a cutoff energy of 8 Ry. Tests were per- 
formed for the electronic structure using a cutoff up to 
15Ry. They did not result in significant changes in the 
character of the wave functions or the formation energies 
of the defects. For the integration in reciprocal space one 
special k-point in the. (2x4) surface Brillouin zone (point- 
group C\h) is used.LJ With this choice of basis set we find 
the theoretical lattice constant to be 5.35 A. nthis value 
is 1.8% smaller than the experimental result,c3 which is 
typical for a III-V semiconductor calculated within the 
LDA (Ref. |24|) and neglecting zero-point vibrations. 

The calculations of the bulk defects are performed in 
a 64-atom cubic supercell. A 96-atom supercell with six 
layers of GaP and a vacuum region corresponding to four 
layers is used for the simulation of defects on and below 
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TABLE I. Formation energies E{ of the gallium surface 
vacancy (Vba) for different cell sizes. All values are given in 
eV and under phosphorus-rich conditions (/xp = 0). The de- 
tails of the calculation are given in Sec. M. 
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the (110) surface. The surface cell in these calculations 
has a (2 x 4) geometry in the [001] and [110] directions 
and the defect is created on one side of the slab. For 
the calculation of defects, structure optimization is per- 
formed to the first two outermost layers for the surface 
defects, to the first three layers for defects in the second 
surface layer, and to nearest and second nearest neigh- 
bors for all other vacancies. In order to quantify the 
errors due to the finite thickness of our slab and surface 
unit cell, we calculated the gallium surface vacancy also 
in larger (3 x 4) and (2 x 5) surface cells as well as for up 
to eight layers in the (2 x 4) cell. These tests resulted in 
a variation of defect formation energies less than 50 meV 
(Table [jj) while the character and eigenvalues of the defect 
states remained almost unchanged. Also the simulation 
of STM images proved to be well converged within these 
variations of cell size and cutoff energy. 

III. SURFACE VACANCIES 
A. The gallium vacancy 

First let us briefly recall the properties of the defect- 
free-GaP(llO) surface, which has been extensively stud- 
iedEj The point-group C\u of this surface consists of a 
single mirror plane in the [110] direction and the relax- 
ation of surface atoms is characterized by a pronounced 
outwards movement of the anion and an inward shift of 
the cation in the surface layer. The result is a character- 
istic buckling of the surface atoms, which is commonly 
described by a tilting angle uj. We find this angle to be 
27° and the height difference between gallium and phos- 
phorus surface atoms to, be.0.56A, in good agreement 
with previous theoreticaOEO and LEED data.ES The re- 
laxation is driven by the lowering of the filled surface 
state located at the surface anions. As a result of the sur- 
face relaxation, the back bonds of the anions are nearly 
perpendicular while the cations show a planar bonding 
configuration indicating that these atoms are p and sp 2 
bonded. The gain in surface energy due to relaxation is 
650 meV per surface unit cell. 

Let us now consider a gallium vacancy Vc a at the sur- 
face: Removing a cation from the surface layer results in 
creating dangling orbitals at the two neighboring phos- 
phorus atoms Pi at the surface and one at the neighbor- 
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FIG. 1. Atomic relaxations for the defects in the first and 
second surface layer. For each defect the top and side view 
of the first two layers are shown. Big atoms belong to the 
surface layer and small ones to the second layer. Crystal- 
lographic directions are indicated as well as the planes taken 
for the side-view pictures (dashed lines crossing the top view), 
(a) Vba in the surface layer; (b) Vp in the surface layer; (c) Vba 
in the second layer; (d) Vp in the second layer. 

ing atom P2 in the second layer of the slab [Fig.|l](a)]. 
The phosphorus atoms Pi in the first layer are twofold 
and the atom P2 in the second layer is threefold coor- 
dinated. The calculated equilibrium positions are sum- 
marized in Table p[ It is interesting to note that only 
the nearest neighbors around the defect significantly re- 
lax. All other atoms remain close to the positions of the 
unperturbed surface with a maximum relaxation smaller 
than 0.1 A. The relaxation of the nearest neighbors is 
characterized by a strong inward relaxation towards the 
vacancy. The atoms Pi move inwards and lift the sur- 
face buckling locally almost completely. We also looked 
in detail for asymmetric distortions by starting from vari- 
ous low-symmetry configurations with equivalent k-point 
sampling. However, for neutral vacancies the energeti- 
cally most stable structure maintains the C\h symmetry 
of the (110) surface. 

Due to the large relaxation of the phosphorus atoms 
towards the center of the vacancy, the distance between 
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TABLE II. Relaxations of the nearest-neighbor atoms of 
the cation (Vba) and anion (Vp) surface vacancy. A x refers 
to the atom of species A at the surface (x = 1) or in the 
second surface layer (x = 2). A[001], A[lT0], and A[110] in- 
dicate the shifts in the directions illustrated in Fig. ^; ||A|| 
is the total displacement, and a the bonding angle of the 
twofold-coordinated surface atoms and of the bonds in the 
plane parallel to the surface for the subsurface atoms. 
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-0.61 
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A[lT0] 


(A) 


0.33 
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A[110] 
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-0.44 


0.27 


-0.29 


0.14 


l|A|| 


(A) 


0.82 


0.56 


0.42 


0.33 
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(deg.) 


90.1 


93.1 


120.3 


98.6 



these atoms is strongly reduced. We find 2.7 A between 
Pi and P2 and 3.1 A for the P1P1 separation, as com- 
pared to 3.8 A for the second nearest-neighbor distance 
in bulk GaP or 2.2 A in black phosphorus t3 This indi- 
cates that the driving force behind the strong relaxation 
is the tendency to form phosphorus-phosphorus bonds. 
For the twofold-coordinated atoms Pi we find a bond 
angle of 90.1°, which is slightly closer to an ideal p con- 
figuration than for the anions at the surface without a 
defect. 

Figure ||(a) summarizes the results for the single- 
particle energies of the gallium surface vacancy. Wc find, 
two defect-related states in the fundamental band gap:c2l 
One that we call la" is located 0.5 eV above the valence 
band, the other called 2a' is close to the upper surface 
band. For the neutral charge state, la" is singly occu- 
pied and 2a' is empty. From analyzing the wave functions 
and the local density of states (LDOS), we find that both 
la" and 2a' consist mainly of p-states centered at the 
three P neighbors. From this analysis we find further a 
resonance la' at about — leV also consisting mainly of 
p-states centered at the three P-neighbors. The lowest- 
lying localized level is determined at about — 8.5eV and 
can be characterized as back bonds of the three nearest- 
neighbor atoms. 

The symmetry and nature of these defect states can, 
be qualitatively understood in terms of a simple modcl:Lil 
Creating a vacancy results in three dangling orbitals at 
the two surface atoms Pi and at the atom P2 in the 
second layer. Rehybridization of these orbitals will give 
rise to three nondegenerated defect states. According 
to the point symmetry of the (110) surface, these de- 
fect states are either symmetric or antisymmetric with 
respect to the mirror plane. This is different from the 
bulk system, where Td symmetry gives rise to ai and ti 
states formed bji_tlie dangling orbitals of four nearest- 
neighbor atoms.EStS Analyzing the character of defect- 
related wave functions, we find the la' state to be sym- 
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FIG. 2. Single-particle energies for surface vacancies, va- 
cancies in the second surface layer, and vacancies in the bulk. 
Top: results for gallium vacancies; bottom: results for phos- 
phorus vacancies. Also indicated are the theoretical values for 
the valence-band maximum VB (E = 0.0 eV), the unoccupied 
surface band SB (1.5 eV), and the conduction-band minimum 
CB (1.7 eV). 

metric while the singly occupied state la" is antisymmet- 
ric. The highest level 2a' is symmetric again. 

To our knowledge this is the first ah initio calculation 
of cation surface vacancies. Allen, Sankey, and Dowo 
calculated the neutral gallium surface vacancy in an em- 
pirical tight-binding scheme: Neglecting atomic relax- 
ations and electronic self-consistency, they found a singly 
occupied level 0.5 eV above the valence-band maximum 
and an empty level at 1.1 eV. These results contrast ours: 
Without including structural optimization, we find two 
nearly degenerate levels above the valence-band maxi- 
mum, but upon atomic relaxation the levels split and 
take the energies shown in Fig.||(a). 

The electronic structure of surface defects can be inves- 
tigated experimentally by scanning tunneling microscopy 
(STM) and scanning tunneling spectroscopy (STS). In 
STM, the charge states can be estimated by examining 
defects on p- and n-type doped crystals and comparing 
the local surface band bending caused by the charged 
defects with that of other charged impurities such as ion- 
ized dopants. EJ In STS, a tunneling spectrum at the de- 
fect site gives directly the position of ionization energies 
in the band gap as has been done for antisite defects at 
the GaAs (110) surfaced From our results, a first esti- 
mate of ionization energies of the defect can be obtained 
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by the positions of the defect levels in the neutral charge 
state. As the gallium surface vacancy has a singly oc- 
cupied state at midgap, we expect at least three charge 
states (+, 0, — ). While we are not aware of STS measure- 
ments on surface vacancies, this result apnears to be in 
contrast to STM measurements by Ebert,E3C3 where the 
gallium surface vacancy was found to be neutral for both 
n- and p-type material. However, it is important to take 
into account that the charge state of a surface defect is 
not determined by the position of the bulk Fermi level 
but by the position of the surface Fermi,-ejuergy. Indeed, 
photoelectron spectroscopy experimentsE3 on GaP(llO) 
indicate that the Fermi level at the surface is pinned 
near midgap, which might stabilize the neutral charge 
state. Further theoretical and experimental studies are 
necessary to clarify this issue. 

B. The phosphorus vacancy 

The phosphorus surface vacancy, Vp, is surrounded by 
three low-coordinated gallium atoms; two in the first 
layer (Gai) and one in the second layer (Ga2). The 
calculated equilibrium geometry of the defect is summa- 
rized in Table|l] and Fig.[|(b). As for the gallium va- 
cancy, the equilibrium structure maintains the C\h sym- 
metry of the defect-free surface and atomic relaxations 
are mainly restricted to nearest-neighbor atoms. Second 
nearest neighbors relax less than 0.1 A. However, the in- 
ward relaxation of the surface atoms Ga! is clearly less 
pronounced than for the gallium vacancy. The origin of 
this is the different configuration of surface cations and 
anions on the defect-free surface: Since surface gallium 
atoms at the relaxed but unperturbed surface are located 
deeper in the surface than phosphorus atoms, the initial 
volume spanned by the nearest-neighbor atoms is much 
smaller for Vp than for Vcja- As a consequence, the relax- 
ation results in an increase of the surface buckling near 
the defect (0.85 A as compared to 0.56 A) and not in a 
decrease as for Vca- 

The distances between the atoms surrounding the de- 
fect are 2.8 A for the GaiGa2 and 3.7 A for GaiGai separ- 
ation. This compares to 3.8 A in bulk GaP and 2.5 A 
in bulk galliumH The angle between the back bonds of 
the surface gallium atoms Gai is 120° as compared to 
122° for the relaxed unperturbed surface. Thus, like the 
cations at the surface without a vacancy, the atoms Gai 
prefer an sp 2 bonding configuration. While this is the 
first ab initio calculation of anion surface vacancies at 
GaP (110), there have been several publications for anion 
vacancies at the (110) surface oi-GaAs: Calculating Va s 
in a (2 x 2) surface cell, Yi et altB determined the height 
difference between surface gallium atoms for the neutral 
vacancy to be 0.2 A and for the negatively charged defect 
to be 0.4 A. Zhang and ZungerEll calculated the inward 
relaxation of the atom Gai towards Va s to be 0.34 A and 
0.29 A in the negative and positive charge state. The 



positively charged vacancy at GaAs CliO) has also been 
calculated by Kim and ChelikowskyEl who give an in- 
ward relaxation of 0.25 A. Keeping in mind that GaP 
has a slightly smaller lattice constant than GaAs, these 
values are in close agreement with our result of 0.3 A. 

The analysis of the electronic structure reveals also 
similarities between cation and anion surface vacancies. 
Again, we find two defect levels in the fundamental band 
gap [Fig. |(d)] . The level la" is located 0.8 eV above the 
valence-band maximum and the level 2a' mixes with the 
unoccupied surface bands close to the conduction-band 
minimum. The occupation numbers are again one and 
zero for these defect states. In the valence band we find 
two localized states at about — 0.8eV and — 6.5eV. Since 
these are resonances in the valence band, they are not 
shown in Fig.[|(d). The levels la" and 2a' consist mainly 
of p-states centered at the gallium neighbors, while the 
lower- lying resonance is mostly s-like in character. As 
for Vcja, the energetically lowest level is symmetric with 
respect to the surface mirror plane, la" is antisymmet- 
ric, and 2a' is symmetric. The electronic structure of the 
anion surface vacancy has been previously calculated by 
Daw and Smithc3e3 using a tight-binding method com- 
bined with Green's- function techniques. Without consid- 
ering atomic relaxations, the#. found the same sequence 
of electronic states for GaAs,c3 namely a symmetric level 
near the valence-band maximum and two levels in the 
band gap, where the lower, antisymmetric level is below 
midgap and the symmetric level is near the bottom of 
the conduction band. In the calculations for GaPjlj the 
authors determined the position of the highest occupied 
level to be at 1.0 eV, which is in good agreement with our 
value of 1.1 eV for the unrelaxed surface vacancy. In con- 
trast to these results, tight-binding calculations by Allen, 
Sankey, and Dowa predicted three deep levels at 1.1, 1.5, 
and 2.2 eV. 

Comparing our results with recent ab initio calcu- 
lations for GaAs (110), we find qualitative agreement 
for the position and character of levels: Kim and Che- 
likowskyE!! obtained for the positively charged vacancy 
the same symmetry of the vacancy states and very sim- 
ilar positions of levels to those find for the neutral 
vacancy at GaP (110). Yi et alx3 determined the half- 
filled level above the center of the band gap while Zhang 
and ZungercJ calculated this level to be 0.41 and 0.5 eV 
above the valence-band maximum for the positive and 
negative charge state. 

Similar to the gallium surface vacancy, we expect also 
for the phosphorus surface vacancy at least two charge- 
transfer levels corresponding to the possible charge states 
(+,0, — ). This rcsult-ls-consistent with STM measure- 
ments by Ebert et a/.,E£lo who determined the phospho- 
rus surface vacancy to be positively charged on p-type 
material and neutral or negatively charged under n-type 
conditions. 
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TABLE III. Relaxations of the nearest-neighbor atoms of 
the cation (Vba) and anion (Vp) vacancy in the third surface 
layer and in the 64-atom bulk supercell. For the vacancies in 
the third surface layer, the relaxations for the nearest neigh- 
bors in the second, third, and fourth layer are given in percent 
of the bulk bond distance (2.32 A). For the bulk vacancies, 
the breathing mode relaxation maintaining the Td symmetry 
is indicated. All values are given in percent of the bulk bond 
distance (2.32 A). 
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IV. SUBSURFACE VACANCIES 

Defects at and near the surface differ from their cor- 
responding bulk defects due to changes in coordination 
number as well as in symmetry. In order to compare 
defects at different positions relative to the surface in a 
systematic way, we have calculated both types of vacan- 
cies in a bulk system as well as in the second and third 
layer of a surface cell. We will first focus on bulk defects 
and vacancies in the third layer. Vacancies in the second 
layer will be discussed separately: They exhibit features 
observed neither for bulk defects nor for surface defects. 



A. Bulk defects 

Bulk vacancies in a zinc-blende semiconductor are 
characterized by a .binding level a\ and a threefold- 
degenerate level t2 .HO'O LJ For vacancies in the neutral 
charge state, the a\ level lies in the valence band and 
the ti states are located above the valence-band max- 
imum for the cation vacancy and near the conduction- 
band minimum for the anionic defect. Table [II sum- 
marizes our results for structural changes and Figs.||(c) 
and ||(f) show the position of levels for the relaxed de- 
fects. Atomic relaxation is mainly governed by an inward 
breathing relaxation conserving the tetrahedral symme- 
try of the lattice. We find the inward relaxation of the 
gallium vacancy to be 3.7% of the bulk bond distance 
(2.32 A). For the unrclaxed vacancy the a\ resonance 
is located — 1.0 eV below the valence-band maximum as 
compared to — 0.75eV calculated by Scheffler et al$3 
with a self-consistent pseudopotential Green's-function 
method. The ti level lies at 0.2 eV, in fair agreement with 
Ref. [?3| (0.15eV). Allowing for atomic relaxation, the a\ 
level shifts to — 1.7eV while the £2 level gains 0.1 eV. 

Calculating the atomic relax atio ns of Vp , we find an in- 
ward relaxation of 6.1% (Table|iTl|). For the ideal atomic 
positions, the a\ resonance is found at about 0.1 eV be- 
low the valence-band maximum. Including atomic relax- 
ations lowers the energy by 0.4 eV. The ti level is found 



1.5 



> 



B° 1 



^ 0.5 



0.0 




0.0 0.2 0.4 0.6 0.8 1.0 



reaction coordinate z 
FIG. 3. Differences in total energy of the vacancies in 

the second layer and the complexes of surface vacancy and 

antisite in the second layer. The reaction coordinate z is the 

lateral position of the nearest-neighbor atom Pi (Gai) of Vba 

(Vp) as denned in Fig. [j](c)(d). At z — 0.0, this atom is 

confined to the plane of the surface anions (cations), and at 

z = 1.0 it is placed in the plane of the cations (anions) in the 

second layer. Solid line, Vgb. and PGa-Vp; dashed line, Vp and 

Gap-VGa- 



within the conduction band for the unrelaxed defect and 
lies below the band edge after relaxation. 

An analysis of our bulk and slab calculations shows 
that these results are almost insensitive to whether the 
vacancies are in a bulk cell or in the third layer of a 
surface cell. Both the positions of the defect levels a\ and 
ti and the formation energies (which will be discussed 
111 Sec.0) agree well. The main differences are slight 
modifications in the atomic geometry. The origin of this 
is the changes in symmetry due to the presence of the 
surface (TableQ). 



B. Defects in the second layer 

Let us now focus on neutral vacancies in the second 
layer. The results for the equilibrium structure are sum- 
marized in Table IV and Figs. |l|(c) and[l](d). As for sur- 
face vacancies, the atomic structure of cation and anion 
vacancies maintains the C\h symmetry of the defect-free 
surface. Analyzing the atomic relaxations, we find signif- 
icant displacements for the nearest-neighbor atom at the 
surface. For both defects this atom moves almost into the 
center of the vacancy. Therefore, vacancies in the second 
layer could also be considered as vacancy-antisite com- 
plexes (PGa-V"p and Gap-Vca)- In contrast to previous 
calculations. £qj bulk vacancies (e.g., Vca and Asc a - Va s 
in GaAs)J£j~t3 we do not find the second-layer subsur- 
face vacancies to be metastable. Figure || shows the en- 
ergy of the defects as a function of the z-coordinate of 
the nearest-neighbor atom Pi (Gai) at the surface. At 
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TABLE IV. Relaxations of the nearest-neighbor and two 
second nearest-neighbor atoms of the cation (Vcj a ) and anion 
(Vp) vacancies in the second surface layer. A x refers to the 
atom of species A in the surface (a; = 1) or second and third 
surface layer (as = 2,3). A[001], A[lT0], and AjllO] indicate 
the shifts in the directions illustrated in Fig. Ill and ||A|| is 
the total displacement. All values are given in A. 
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z = 0.0 this atom is confined to the plane of the surface 
anions (cations), and at z = 1.0 it is placed in the plane 
of the cations (anions) in the second layer. In the equi- 
librium structure, we find the total displacement of the 
atom Pi next to Vq & to be 1.8 A, and that of the atom 
Gai next to Vp is 1.0 A. These values can be compared 
to 2.3 A for the bond length between atoms in the first 
and second layer on the relaxed unperturbed surface. In 
contrast to the surface vacancies, we find a substantial 
displacement also for the second nearest-neighbor atoms 
in the surface layer. The two nearest-neighbor atoms in 
the second and the one in the third layer of the slab un- 
dergo almost no relaxation (Table [V ) . 

The gallium vacancy has two levels within the fun- 
damental band gap [Fig.||(b)]: one doubly occupied at 
0.6 eV and one singly occupied at about 1.2 eV. These 
states are antisymmetric and symmetric with respect to 
the surface mirror. For the phosphorus vacancy, we find a 
symmetric and singly occupied defect state at the middle 
of the band gap and an empty and antisymmetric state 
below the conduction-band minimum [Fig.^(e)]. As for 
surface vacancies, we expect the cation and the anion va- 
cancy in the second layer to show an amphoteric behav- 
ior: While under p-type conditions they should act as 
a donor, n-type conditions might stabilize a negatively 
charged state. Thus cation and anion vacancies in the 
first and second layer are possible candidates for pinning 
the surface Fermi level. 

We conclude that the equilibrium structure of sub- 
surface vacancies is mainly determined by the coordina- 
tion of the nearest-neighbor atoms: Twofold-coordinated 
atoms in the surface layer can alter their position signifi- 



TABLE V. Formation energies for surface and subsurface 
vacancies under phosphorus-rich conditions. The table shows 
the results for unrelaxed and relaxed defects per broken bond 
in order to compare the threefold-coordinated surface vacan- 
cies to the fourfold-coordinated subsurface and bulk defects. 
All values are given in eV. 
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cantly without breaking the remaining bonds. However, 
threefold-coordinated atoms as completely surrounding a 
vacancy in the third or a deeper layer are strongly bound 
to their positions in the unperturbed lattice. Thus the 
ability of atoms in the second layer to relax is compara- 
ble to that of bulk atoms. Consequently, vacancies with 
only threefold neighbors (i.e., vacancies below the second 
layer) show very similar properties to bulk defects, while 
vacancies that have twofold-coordinated neighbors (i.e., 
vacancies at the surface and in the second layer) show a 
very different geometry and electronic structure. 

V. DEFECT FORMATION ENERGIES 

In the preceeding sections we have focused on the 
atomic geometry and the electronic structure of the de- 
fects, and how these properties are modified if a bulk de- 
fect is brought to the surface. We will now discuss how 
the energy required to create a vacancy in thermody- 
namic equilibrium (i.e., the formation energy) is affected 
by the surface. The formation energy of a defect is given 
by 

Et = £ t d f ct - AWca - N Pl i P - qE Fcrmi . (1) 

A^Qa and Np are the number of cations and anions with 
chemical potentials /^Ga and /ip, and q is the number of 
excess electrons in exchange with the electron reservoir 
^Fermi- Since we are focusing on neutral defects (q = 0), 
the formation energies are independent of the Fermi level. 
The Ga and P chemical potentials are related to each 
other by the condition of chemical equilibrium with the 
bulk GaP crystal: 

/i G aP = MGa + MP ■ (2) 

The formation energy thus can be written a E 

E{ = ^defect _ ^ {Np _ NGa) _ AT GaM bulk _ (3) 
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Moreover, we can approximate the upper limit of the 
chemical potential of each species by their corresponding 
bulk phases (Hc^ k , /i P ulk ). With the definition of the 
heat of formation of the crystal, 



A tt ..bulk ..bulk hi 

Ai/f - ^ GaP - n Ga - /i P 



bulk 



bulk 



(4) 



we get a range in which the chemical potential can be 
varied from gallium-rich to phosphorus-rich conditions: 



bulk 



AH{ < hp < up 



bulk 



(5) 




and black 
to be the equi- 



Assuming metallic orthorhombic galliU| 
phosphorus in an orthorhombic unit eel 
librium structures of the condensed elemental phases, we 
get AH{ = 1.13 eV. This value-is slightly higher than the 
experimental result (1.08 eV) £3 

Figure |] summarizes our results for the different re- 
laxed defects. Over a wide range of the chemical poten- 
tial, the anion surface vacancy is the energetically most 
favored defect. Only under very phosphorus-rich condi- 
tions is the cation surface vacancy preferred. Vacancies 
in the second and third layer show a significantly higher 
formation energy. Creating a gallium vacancy in the sec- 
ond layer takes 0.25 eV more than at the surface while 
the other subsurface vacancies are even less favored. 

Table lists the formation energies with and without 
considering atomic relaxations. Note that all energies in 
Table are given per broken bond and not per defect. 
This allows us to compare directly bulk and surface va- 
cancies that have different coordinations. For simplicity, 
all formation energies are given under phosphorus-rich 
conditions. According to Eq. the following discus- 
sion remains unaffected by considering a different chemi- 
cal environment. Let us first discuss the case without any 
atomic relaxation: The most striking feature is the fact 
that the energy per broken bond is almost unchanged 
with respect to-the position of the vacancy (surface, sub- 
surface, bulk).E3 The deviations are smaller than 0.1 eV, 
i.e., close to the estimated error of our method. The 
largest difference (70meV) is found for Vba between the 
first and second surface layer. 

However, if we consider atomic relaxation, the picture 
changes qualitatively: For vacancies in the first and sec- 
ond layer, we find a strong relaxation. Due to this re- 
laxation the formation energy is largely reduced (up to 
0.5 eV, i.e., the energy gained by atomic relaxation is one 
order of magnitude larger than for bulk vacancies). It is 
further interesting to note that the relaxation energy for 
Vba is about twice as large as for Vp. 

For bulk vacancies and vacancies in the third surface 
layer, the lowering of the formation energy due to atomic 
relaxation is small and almost in depen dent of the type of 
defect. It has been shown in Sec. IV A that these defects 



are very similar with respect to defect levels and atomic 
equilibrium positions. The presence of the surface also 
induces differences in formation energy that are less than 
40meV indicating the bulklike nature of the subsurface 
vacancies. 
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FIG. 4. Formation energies for relaxed surface and sub- 
surface vacancies as a function of the difference of the chemi- 
cal potential of phosphorus and bulk phosphorus. The limits 
for Ga-rich (hp — Hp u = — 1.13 eV) and P-rich conditions 
(hp - /ip ulk = 0.00 eV) are defined in the text. Note that the 
defects labeled as Vba (second layer) and Vp (second layer) are 
better described as Gap - Vp and PGa - Vba complexes. 



VI. SIMULATION OF STM IMAGES 
A. Method 

For simulating S TM . images we employ the Tersoff- 
Hamann approach.! 5 Tr 8 ! In this approach the tunneling 
current is proportional to the local density of states at 
the tip position, p s (r tlp ,e), integrated over the interval 
between the Fermi surface and the applied bias U: 



I(r tip ,U) oc 



(6) 



A positive bias results in imaging the empty states 
in the specified energy range at position r tlp and nega- 
tive bias in imaging the filled states. By mapping sur- 
faces of constant I(r tlp ,U) an isocurrent picture can be 
simulated. Alternatively the integration might be evalu- 
ated for fixed tip height. Calculations of various recon- 
structed semiconductor surfacesKTjM.as well as surfaces 
with impurities and adsorbatescJoHa have proven this 
approach to be very useful for the interpretation of ex- 
perimental results. 



B. Results 

Figures ||(a) and ||(b) display our results for the oc- 
cupied states of the gallium and phosphorus surface va- 
cancy. The plots were taken at a distance of 3 A from 
the outermost anion layer and include all occupied defect 
and surface states as well as surface resonances down to 
1.4 eV below the valence-band maximum. The pictures 
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FIG. 5. STM simulations of occupied states for sur- 
face and subsurface vacancies, (a) Gallium surface vacancy; 
(b) Phosphorus surface vacancy; (c) Vgb, in the second layer; 
(d) Vp in the second layer; (e) Vg& in the third layer; (f) Vp in 
the third layer. The electron density in the plots varies from 



0.4xl(T 7 A to4.0xl(T 
is estimated to be 3 A. 



and the tip surface separation 



are dominated by the filled dangling bonds above the an- 
ions. Thus, for the defect-free surface the STM picture 
shows the periodicity of the surface lattice and the max- 
ima of the charge density are close to the positions of the 
surface anions. 

If we focus on the cation vacancy [Fig.||(a)] we find a 
localized perturbation in the lattice: the dangling bonds 
of the two neighboring atoms along the [110] direction 
are laterally displaced according to the relaxation of the 
atoms and clearly appear less intense than the ones far- 
ther away from the defect. This signature of_th£_defect 
has also been identified in STM experiments. How- 
ever, from the experimental data it is not clear whether 
this effect is due to atomic relaxation or due to changes in 
the electronic structure near the vacancy. We therefore 
compared the LDOS above the surface dangling bonds 
of the neighboring atoms with that of dangling bonds 
farther away from the defect. The results are shown in 
Fig. ^| together with an inset defining the lateral positions 
where the tip height has been varied. In the investigated 
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FIG. 6. The local density p s as a function of the tip 
height above the surface. The inset shows the same part of the 
surface as Fig.|^(a) (gallium surface vacancy). Squares corre- 
spond to the position of the nearest-neighbor surface dangling 
bond and circles to that of the next anion in the direction of 
the zigzag row. The decay constant deduced from fitting the 
data is 3.1 A" 1 for both graphs and the height difference for 
constant density is 0.42 A. As in Fig.|j| p a includes all occu- 
pied defect and surface states as well as surface resonances 
down to 1.4 eV below the valence-band maximum. The zero 
of the abscissa corresponds to the position of the outermost 
anions at the relaxed ideal surface. 



interval above the surface we find an exponential decay 
of the density into the vacuum. The decay constants 
at the two maxima are virtually identical (3.1 A -1 ). We 
therefore expect that in constant current STM the height 
difference between the two maxima is not a function of 
the tunneling current. From Fig. ^ we find this height dif- 
ference to be 0.42 A, which is very close to the difference 
in the atomic structure (0.44 A, Table||). We note that 
the good agreement between the height corrugation in 
the electronic and atomic structure is related to the spe- 
cific structure of the (110) surface and that of the defect: 
Both rehybridization and charge transfer for this system 
are small, and our STM simulation mainly shows the 
filled surface dangling bonds, which are not significantly 
changed by the presence of the vacancy. This is true 
even for nearest-neighbor atoms. We expect larger de- 
viations for charged surface defects or for more complex 
systems where defect states are less localized or where 
rehybridization of atoms near the defect is more promi- 
nent. 

The STM simulation of the anion surface vacancy is 
characterized by a depletion of the density of states due 
to the missing dangling bond [Fig.||(b)]. The perturba- 
tion is very localized above the defect site and already 
the dangling bonds located at the second nearest neigh- 
bors are not influenced by the defect. This .-agrees well 
with experimental data for the neutral defecuIS and with 
the short-range character of the defect's perturbation of 
the surface as discussed in Sec. HI. 
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Simulating the vacancies in the second layer, we get the 
results shown in Figs.[|(c) and ||(d). Due to the large re- 
laxation of the anion in the surface layer, the gallium va- 
cancy displays features similar to the surface phosphorus 
vacancy: The phosphorus atom above the vacancy has no 
states within the simulated energy range at that location, 
which can be explained by the large vertical displacement 
of this atom. For the phosphorus vacancy the dangling 
bonds of the second nearest neighbors are slightly dis- 
placed and depressed. As for the cation surface vacancy, 
we analyzed this behavior as a function of the tip height. 
We find the height difference of the tip to be 0.22 A as 
compared to 0.27 A structural relaxation. Thus, also the 
signature of this vacancy is mainly governed by the re- 
laxation of the surface atoms. We note that the STM 
simulations of the gallium surface vacancy [Fig.^(a)l and 
of the phosphorus vacancy in the second layer [Fig.g(d)] 
show a remarkable resemblance. The same holds for the 
anion surface vacancy [Fig.^](b)] and the cation vacancy 
in the second layer [Fig.||(c)]. Therefore, we expect these 
defects to be hardly discernible in STM measurements of 
the occupied states. However, for the gallium vacancy at 
the surface and Vp in the second layer, the highest oc- 
cupied defect state has a different symmetry [Figs.^(a) 
and ||(e)]. Thus, both defects should be distinguishable 
in experiment. Also Vp at the surface and Vba in the sec- 
ond layer should be distinguishable in experiments such 
as scanning tunneling spectroscopy because they exhibit 
very different positions of the highest occupied defect 
state [Figs. |(b) and|(d)]. 

The simulation of a gallium vacancy in the third layer 
is virtually identical to that of a defect-free surface 
[Fig.|5|(e)]. For the phosphorus vacancy, a slight depres- 
sion above the missing atom is found [Fig.^(f)], which, 
however, might be too weak to be resolved in experi- 
ment. We therefore expect that neutral vacancies be- 
low the second surface layer cannot be detected by STM 
measurements. For charged defects the local surface band 
bending induced by the impurity might change this result 
quahtativelyfm 

VII. SUMMARY 



surface layer on, the defects show essentially the same 
properties as bulk vacancies, namely a breathing- mode- 
like relaxation of the nearest-neighbor atoms and a triply 
degenerate acceptor (Vba) and donor level (Vp), respec- 
tively. Vacancies in the second surface layer are found 
to be unstable against the formation of vacancy anti- 
site complexes where the nearest-neighbor surface atom 
moves towards the vacancy site. For Vca, the highest 
occupied level lies near the empty surface band while for 
Vp it is found at a midgap position. 

As can be expected from their lower coordination num- 
ber, surface vacancies are energetically preferred as com- 
pared to subsurface and bulk vacancies. Our results 
show that, especially at and near the surface, the forma- 
tion energy of defects is significantly affected by rather 
large atomic relaxations. For all vacancies with nearest- 
neighbor atoms in the surface layer we find relaxation en- 
ergies that are at least six times higher than for defects 
with four threefold-coordinated nearest-neighbor atoms. 

Our calculated STM images of the surface vacancies 
agree well with available experimental data. The com- 
parison of the STM simulation with the atomic struc- 
ture shows that neutral vacancies at this surface exhibit 
a close relation between atomic relaxation and height 
corrugation in the STM picture. To our knowledge, 
unlike charged subsurface defects, which have been fre- 
quently observed using STM, uncharged subsurface va- 
cancies have not been reported by experimental groups so 
far. This is consistent with our results for vacancies from 
the third layer on, which exhibit virtually no interaction 
with the surface. However, vacancies in the second layer 
look like surface vacancies in our simulation. In experi- 
ment they might be discerned from surface vacancies by 
spectroscopic methods. 
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Using first-principles total-energy calculations we in- 
vestigated neutral surface and subsurface vacancies on 
both sublattices at GaP (110). Atomic relaxations, elec- 
tronic structure, as well as formation energies were an- 
alyzed, and STM pictures of the various vacancies were 
simulated. 

For both the gallium and the phosphorus surface va- 
cancy, we find a pronounced inward relaxation of near- 
est neighbors into the defect. The electronic structure 
is characterized by a singly occupied deep level for Vc a 
and Vp. Thus we expect them to be amphoteric, indi- 
cating that surface vacancies are electrically active and 
potentially a source of compensation. From the third 
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